A Novel Mode Switching Scheme Utilizing Random 
Beamforming for Opportunistic Energy Harvesting 

Hyungsik Ju and Rui Zhang 



Abstract 

Since radio signals carry both energy and information simultaneously, a unified study on simultaneous wireless 
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information and power transfer (SWIPT) has recently attracted a great deal of attention for providing wireless 



energy to communication devices. This paper proposes a novel receiver mode switching scheme for SWIPT and 
the enabling transmitter design for a point-to-point multiple-input single-output (MISO) channel when the channel 
state information (CSI) is not available at the transmitter. The proposed scheme employs random beamforming at 



the transmitter to generate artificial channel fading to enable opportunistic energy harvesting at the receiver when 
the channel power exceeds a certain threshold. For the proposed scheme, we investigate the achievable average 
^ ■ information rate, average harvested energy, and energy outage probability, as well as their various trade-offs in 

fading channels. Compared to a reference scheme of periodic receiver mode switching without random transmit 
beamforming, the proposed scheme is shown to be able to achieve better rate-energy trade-offs when the harvested 
energy target is sufficiently large. Particularly, it is revealed that employing one single random beam for the proposed 
scheme is asymptotically optimal as the transmit power increases to infinity, and also performs the best with finite 
transmit power for high energy harvesting requirements of practical interests, thus leading to an appealing low- 



complexity implementation. Finally, we compare the rate-energy performance of the proposed scheme with selected 



random beam designs. 



Index Terms 



Simultaneous wireless information and power transfer (SWIPT), wireless power, mode switching, multi-antenna 
system, random beamforming, rate-energy trade-off, energy outage. 

I. Introduction 

Conventionally, fixed energy supplies (e.g. batteries) are employed in energy-constrained wireless networks, 
such as sensor networks. The lifetime of the network is typically limited, and is thus one of the most important 
considerations for designing the network. To prolong the network's operation time, energy harvesting has recently 
attracted a great deal of attention since it enables scavenging energy from the environment and potentially provides 
unlimited power supplies for wireless networks. 

Among other commonly used energy sources (e.g. solar and wind), radio signals radiated by ambient transmitters 

have drawn an upsurge of interest as a viable new source for wireless energy harvesting. Harvesting energy 
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Fig. 1. A MISO wireless system for SWIPT via receiver mode switching. 

from radio signals has already been successfully implemented in applications such as passive radio-frequency 
identification (RFID) systems and body sensor networks (BSNs) for medical implants. More interestingly, wireless 
energy harvesting opens an avenue for the joint investigation of simultaneous wireless information and power 
transfer (SWIPT) since radio signals carry energy and information at the same time. SWIPT has recently been 
investigated for various wireless channels, e.g., the point-to-point additive white Gaussian noise (AWGN) channel 
IH, O, the frequency selective and/or fading channels |3), 0, the multi-antenna channels (5J, Q, Q, and the 
multiple-access and relay channels JS). 

To achieve optimal wireless energy transfer (WET) and wireless information transfer (WIT) simultaneously, 
one key challenge is to develop receiver architectures to enable information decoding (ID) and energy harvesting 
(EH) from the same received signal at the same time [2]. Practically, two suboptimal receiver designs for SWIPT 
have been proposed in [5] based on the principle of orthogonalizing ID and EH, namely power splitting and time 
switching. The power splitting scheme splits the received signal into two streams of different power for ID and 
EH separately, while the time switching scheme switches the receiver between an ID mode and an EH mode from 
time to time. The optimal switching rules between ID versus EH modes for a point-to-point single-antenna fading 
channel subject to the co-channel interference have been derived in H to maximize the average harvested energy 
given an information-rate target. It was shown in [)4| that the time-fluctuation or fading of wireless channels is 
indeed beneficial for receiver mode-switching (time-switching) based SWIPT systems, where an "opportunistic" 
energy harvesting scheme is proved to be optimal, i.e., the receiver should switch to the EH mode when the channel 
power is larger than a certain threshold, and to the ID mode otherwise. 

In this paper, we further investigate the receiver mode-switching based SWIPT in a point-to-point multiple- 
input single-output (MISO) system as shown in Fig. [TJ It is assumed that the MISO channel from the transmitter 
to the receiver is quasi-static flat-fading, i.e., the channel remains constant during each block transmission, but 
varies from one block to another. It is also assumed that the MISO channel is perfectly known at the receiver in 
each block but unknown at the transmitter. Furthermore, we assume that the receiver has no given power supply 
and needs to replenish energy from the received signals via WET based on time-switching, i.e., the receiver can 



either decode information or harvest energy from the received signals at any time, but not both. In conventional 
wireless communication systems for WIT only, multi-antenna techniques have been widely applied to achieve spatial 
multiplexing (9] and/or spatial diversity |[T0l , iTTTI . Furthermore, the use of multiple transmit antennas has been 
considered in wireless MISO broadcast channels to enable opportunistic communication by employing random 
transmit beamforming [12|, ||T3l , which generates artificial time-variations in quasi-static channels. In this paper, 
we apply the random beamforming technique in MISO point-to-point systems for a new usage, i.e., to enable 
opportunistic energy harvesting at the time-switching receiver for more efficient SWIPT implementation. 
The main contributions of this paper are summarized as follows: 

• We propose a novel receiver mode switching scheme for SWIPT based on random transmit beamforming. 
The proposed scheme employs random beamforming at the multi-antenna transmitter to generate artificial 
channel fading in each transmission block during which the physical MISO channel is assumed to be constant. 
This is to most efficiently utilize the time-switching receiver for optimizing the simultaneous WET and WIT 
performance. For the proposed scheme, we characterize the performance trade-offs between WET and WIT by 
investigating the achievable rate and harvested energy in each transmission block given a fixed MISO channel, 
assuming Gaussian distributed random beams. 

• We extend our analysis for the constant MISO channel to quasi-static fading MISO channels. We investigate 
the achievable average information rate and average harvested energy over different transmission blocks, and 
characterize their asymptotic trade-offs when the transmit power goes to infinity. It is shown that employing 
one single random beam for the proposed scheme achieves the best rate-energy (R-E) trade-off in fading 
MISO channels asymptotically. Furthermore, we compare the R-E performance of our proposed scheme with 
that of a reference scheme with periodic switching between ID and EH modes, but without random transmit 
beamforming applied. 

• When the receiver consumes significant amount of energy at each block and/or the capacity of its energy 
storage device is limited, the receiver may suffer from energy shortage unless the amount of harvested energy 
in each block is larger than a certain requirement. We thus study the "energy outage probability" of the 
proposed scheme in fading MISO channels, which is compared to the periodic switching in both asymptotic 
and finite transmit power regimes. 

• In practice, transmit power is preferable to be constant for the optimal operation efficiency of transmitter 
amplifier. However, the use of Gaussian distributed random beams for the proposed scheme can cause large 
transmit power fluctuations. We thus propose alternative random beam designs with more stable transmit power, 
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for which the R-E performance is characterized and compared with the case of Gaussian random beam. 

The rest of this paper is organized as follows. Section |TT] introduces the proposed scheme, and derives the 
harvested energy and achievable information rate for a block of given MISO channel. Section [III] investigates the 
R-E performance for the proposed scheme in fading MISO channels. Section [TV] compares the performance of the 
proposed scheme with different random beam designs. Finally, Section [V] concludes the paper. 

Notations: In this paper, matrices and vectors are denoted by bold-face upper-case letters and lower-case letters, 
respectively. Ijv denotes an N x N identity matrix and represents a matrix with all zero entries. The distribution 
of a circularly symmetric complex Gaussian (CSCG) random vector with mean vector fi and covariance matrix £ 
is denoted by £A/"(/x, X), and ~ stands for "distributed as". C axb and M, axb denote the spaces of a x b matrices 
with complex and real entries, respectively. ||z|| denotes the Euclidean norm of a complex vector z. E [•] represents 
the statistical expectation. 

II. System Model 

As shown in Fig. \T\ this paper considers a point-to-point SWIPT system consisting of one transmitter (Tx) and 
one receiver (Rx). Tx is equipped with N t > 1 antennas and Rx is equipped with one single antenna. It is assumed 
that the MISO channel from Tx to Rx follows quasi-static flat-fading, where the channel remains constant during 
each block transmission time, denoted by T, but varies from one block to another. It is further assumed that the 
channel is perfectly known at Rx, but is not known at Tx. 

The transmitted signal at the nth symbol interval in one transmission block of interest is denoted by x (n) G 
C NtXl . The covariance matrix of the transmitted signal is thus given by S x = E[x (n) x H (n)] = ;^-I/v t > where 
P denotes the total transmit power, which is assumed to be equally allocated among N t transmit antennas. The 
received signal at Rx is then expressed as 

y (n) = h T x (n) + z (n) , (1) 

where y (n) and z (n) denote the received signal and noise at Rx, respectively; and it is assumed that z in) ~ 
CM (0,1) and is independent over n. The MISO channel from Tx to Rx is denoted by h G C NtXl , which is 
constant during each transmission block. 

In each block, Tx aims at achieving SWIPT to Rx. It is assumed that Rx is equipped with a rechargeable battery 
to store the energy harvested from the received signal, which is used to provide power to its operating circuits. 
Specifically, Rx harvests energy from the received signals when it is in the EH mode, while it decodes information 
in the ID mode. We assume that Rx switches between ID mode and EH mode as in H, (H since it is difficult yet 
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Fig. 2. Transmitter and receiver structures for periodic switching (PS). 
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to use the received signal for both ID and EH at the same time due to practical circuit limitations ||2]. As in [@], 

ID mode and EH mode are represented by defining an indicator function as 

{1, ID mode is active 
(2) 
0, EH mode is active. 

We consider two time switching schemes, namely "periodic switching (PS)" and "threshold switching (TS)" as 
elaborated next. 

A. Reference Scheme: Periodic Switching 

As shown in Fig. El with PS Rx sets p = 1 during the first tT amount of time in each transmission block, with 
< r < 1, and p = for the remaining block duration (1 — t)t13 For given h and r, the amount of harvested 

energy (normalized by T) in a transmission block can be derived using S x as 

£(P) (H,t) = (1- t)(E [||h T x(n)|| 2 
= (1 - r) (PH, 

1 1 1 2 

where H = ||h|| is the average channel power per transmitting antenna, and < ( < 1 is a constant reflecting 
the loss in the energy transducer when the harvested energy is converted to electrical energy to be stored. In d3), it 
has been assumed that the energy harvested due to the receiver noise is negligible and thus is ignored. It is further 
assumed that ( = 1 in the sequel for notational brevity. 

The transmitter structure for PS is also shown in Fig. |2] Note that with PS Rx can adjust r based on its energy 
and rate requirements, as well as the channel condition. However, instead of requiring Rx to feed back the values 
of r to Tx to maintain synchronized information transmission, which incurs additional overhead, Tx can employ 
the rateless code Ifl4l for channel coding to realize a seamless information transmission to Rx regardless of to The 

'ideally, with a given time allocation r, setting p = 1 or at the beginning of each block will not change the system performance; 
however, setting p = 1 initially may be practically more favorable for Rx to implement block-wise time synchronization. 

2 This is especially useful when the setup of this paper is extended to the case of multiple receivers in a wireless information-multicast 
channel, where receivers can set different values of r for decoding common information sent by the transmitter, based on their channel 
conditions and energy requirements. 
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bit stream to be transmitted during a transmission block is thus first encoded by a rateless code. Space-time (ST) 
code is then applied to modulate the output bits from the rateless-code encoder, and the modulated symbols are 
transmitted by N t antennas. We consider a ST code of length L, denoted by matrix X( p ) G C LxNt . It is assumed 
that X( p ) is a capacity-achieving ST code[] 

Tx transmits a sequence of X( p )'s in each transmission block. Considering X( p ) with L consecutive transmitted 
symbols from each antenna, (fl} is modified as 

y = X< p >h + z, (4) 

where y G C Lxl and z € C Lxl denote the received signal vector and noise vector, respectively, and z ~ CM (0, IfJ. 
Since X( p ) is assumed to be a capacity-achieving ST code, the achievable rate of the channel in (H} can be shown 
equivalent to that of a MISO channel h with input covariance matrix S x = ^-IjVf Assume that the number of ST 
coded blocks transmitted in each block is sufficiently large such that tT is approximately an integer number of the 
ST block durations for any value of r. For given h and r, the information rate for PS can thus be expressed as 

R<V(H,r)=Tlog 2 (l + PH). (5) 

Note that i?( p ) (H, t) is achievable when N t < 2, but is in general an upper bound on the achievable rate when 
N t > 2 for given h and r. 

B. Proposed Scheme: Threshold Switching 

As shown in Fig. [3l the TS scheme is designed to take advantage of the received signal power fluctuations 
induced by transmit random beamforming within each transmission block for EH/ID mode switching, even with 
a constant MISO channel h. For this purpose, each transmission block is further divided into K sub-blocks and 
artificial channel fading over different sub-blocks is generated by employing multi-antenna random beamforming 
at Tx. 

Furthermore, at the fcth sub-block, k = 1, . . . ,K, Rx determines whether to switch to ID mode or EH mode 
based on A{k), which denotes the channel power at the fcth sub-block (to be specified later). According to |@], 
in the presence of received channel power fluctuations, the optimal mode switching rule that achieves the optimal 
trade-off between the maximum harvested energy and information rate in a transmission block is given by 

, , f 1, if A (k) < A 
p(k) = { ' { } ~ (6) 

[ 0, otherwise, 

3 Alamouti code |10| is known as the capacity-achieving ST code when Nt = 2. For N t > 2, capacity-achieving ST code has not yet 
been found in general. In this paper, however, capacity-achieving ST code is assumed even when N t > 2 to provide a performance upper 
bound for the system under consideration. 
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Fig. 3. Transmitter and receiver structures for threshold switching (TS). 

where A > is a pre-designed threshold on the channel power. It is noted that choosing EH or ID mode at the 
kth sub-block is determined by the channel power A(k), and the threshold on the corresponding received signal 
power thus increases linearly with P, which is given by PA. 

Artificial channel fading over sub-blocks is generated at Tx by using N random beams simultaneously, 1 < N < 
N t . Denote the nth random beam (RB) at the kth sub-block as cf> n (k) G C NtXl , where 4> n (k) ~ CM(0, ^IjvJ 
assuming Gaussian RBs[j and K[4> n (k) 4>mU)] = if k j and/or n ^ m. Then it follows that A (k) = 



TvH a 



! , where a(/c) = $> T (k)h 6 



->Nxl 



is the equivalent MISO channel at the kth sub-block generated by 



$ (k) = [^(k) (f> 2 (k) ■ ■ ■ 4> N (k)}. It can be easily verified that a(k) ~ CM (0, H1 N ) for a given H, and A (k) 
is thus a chi-square random variable with 2N degrees-of-freedom. With N RBs and conditioned on a given MISO 
channel realization h with H = h, the probability density function (PDF) of A := A(k), \/k, and the cumulative 
distribution function (CDF) of A are given, respectively, as [15] 

1 



{h/N) N r (n) 



_ a N-l e -(N/h)a 



F ( "ha\h) . 

I T(N) 



(7) 
(8) 



where F (x) = J °° t x ~ 1 e~ t dt and T (a, x) = t a ~ 1 e~ t dt represent the Gamma function and incomplete Gamma 
function, respectively. From d7J and d8j, it can be inferred that the artificial channel fading induced by random 
transmit beamforming is more substantial over time with smaller N. 

Similarly to PS, Tx with TS generates ST modulated symbols to be transmitted through the artificial fading MISO 
channels. However, they differ in two aspects. First, Tx with TS should employ the erasure code [16] instead of the 
rateless code in the case of PS. This is because the set of sub-blocks used for ID according to © are in general 



Alternative RB designs will be studied later in Section IIVI 



randomly distributed within a transmission block with A > 0, and thus the resulting channel from Tx to Rx in ID 
mode can be modeled by an erasure channel for TS. Second, the ST code is applied over N RBs with TS instead 
of N t antennas with PS. This is because the use of N RBs transforms the JVt x 1 constant MISO channel h into 
an N x 1 fading MISO channel specified by a(fc)'s in each transmission block. For each of the K sub-blocks in 
TS, we consider the use of a ST code of length L denoted by matrix X^ T ) G C LxN . For convenience, we express 
X( T ) = [x^ • • • x^] T , where x| T ^ E C JVxl , 1 < I < L, denotes the Ith transmitted signal vector in each 
ST coded block. The covariance matrix for x z is given by SW = E[x, (x , ) ] = -j^Ijv, V/, to be consistent 
with S x = jjf-IiVf Similar to X( p ) in the case of PS, X^ T ) is assumed to be a capacity-achieving ST code for an 
equivalent MISO channel with N transmitting antennas. 

The received signal at each sub-block is used for either energy harvesting or information decoding according to 
((6]). For the kth sub-block, the received signal can be expressed by modifying (Q]) as 

y = X( T )$ T (£;)h + z 

(9) 

= X( T )a(A:) + z, 

where y G C Lxl and z G C Lxl denote the received signal and noise vectors, respectively, with z ~ CA/"(0,Il). 
When p (k) = 0, the amount of harvested energy (normalized by sub-block duration T/K) at the Mi sub-block is 



derived using { as 



1j 



PA(k). (10) 



Furthermore, by assuming a capacity-achieving ST code, the achievable rate with TS at the /cth sub-block when 
p(k) = 1 can be expressed as 

R^(k) = log 2 (l + PA(k)). (11) 

The amount of harvested energy in a transmission block is the sum of the energy harvested from all sub-blocks 
in the EH mode. Assuming that K — > oo, the harvested total energy (normalized by T) in a transmission block for 
given N RBs, threshold A, and the realization of the MISO channel h with H = h can be obtained from (fTOl as 

£ (T) (h, N,A) = ± lim £ (1 - P (A;)) T X E ™ {k) 

k=l 

= E[(1 -p(k))PA(k)] 

/•oo 

= / Pa A N H (a\h)da (12) 
J A 

t ( n + 1j nA\ 

= Ph^ (13) 

r(jv + i) ' UJ 
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where (Ql) can be obtained by applying © and El 3.351-2] to (O. From £l3]>, it is observed that E^ (h, N, A) 
increases linearly with P. Thus, for a given P, the energy harvesting performance in a transmission block can be 
measured by the normalized harvested energy, £( T ) (h,N,A) jP. Fig. 1 shows E^ (h,N,A) jP versus different 
values of h for P = 20dB when N = 1, 2 and A = 0.1, 0.2, 0.5. It is observed that E^ (h,N,A) jP decreases 
with increasing A when N and h are both fixed, which is in accordance with dT3l . Moreover, when N and ^4 are 
both fixed, E^ (h,N,A) jP is observed to increase monotonically with h. This is because i 7 ^^ {A\h) in ([U) 
decreases with increasing h, and thus 1 — (A \ h), which is the percentage of the received sub-blocks allocated 
to EH mode in each block, increases. Thus, the amount of harvested energy in each block increases with h thanks 
to the increased number of sub-blocks assigned to EH mode, as well as the increased average channel power h, as 
can be inferred from (Q~3}. 

Furthermore, when h and A are both fixed, E( T > (h, N,A) / P is observed to decrease with increasing N when 
h is small, but increase with N when h is sufficiently large. This is because the artificial channel fading is more 
substantial when smaller number of RBs, N, is used, although the same average channel power is given as h. 
Given 1 < N < Nt, it can be shown that F^ (A\h) in © increases with N when h is small, and thus larger 
energy is harvested with smaller number of RBs. In contrast, it can also be shown that F^^ (A\h) decreases with 
increasing N when h is larger than a certain threshold, and thus more energy is harvested with larger number 
of RBs. Similarly, we can verify that E^ (h, N, A) /P increases with N when A is small, but decreases with 
increasing N when A is sufficiently large. 
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Next, the achievable rate in a block for given N, A, and h can be derived from © and (fTTT i as 

(h, N, A) = E [p (k) log 2 (1 + PA (&))] 
/•j4 

= J log 2 (1 + Pa) (a \h)da. (14) 

With f^zj (a \h) given in ([7]), it is in general difficult to obtain a unified closed-form expression of (fl4l for arbitrary 
values of N. However, it is possible to derive closed-form expressions for (fT4l for some special values of N. For 
example, (h, 1, A) and R^ (h, 2, A) for iV = 1 and 2, respectively, are derived in closed-form in Appendix 
El Fig. shows i?( T ) (/i, iV, A) versus different values of h f or P = 20dB when iV = 1, 2 and A = 0.1, 0.2, 0.5. It 
is observed that i?( T ) (/i, N, A) increases with A when N and /i are both fixed, which is in accordance with (fl4l . 
Moreover, by the opposite argument of the explanation for Fig. 01 when h and A are both fixed, R^ (h, N, A) 
is observed to increase with iV when h is small or A is large, but decrease with increasing N when h or A is 
sufficiently large/small. 

However, different from {h,N,A^ /P in Fig. @] which is a monotonically increasing function of h, it is 
observed in Fig. |5] that R^ (h, N, A) in general first increases with h, and then decreases with increasing h for 
given N and A. The reason is as follows. When h — > 0, from (fl4l . we have i?( T ) (/i, TV, A) -»• log 2 (1 + Ph); thus, 
i?( T ) (h,N, A) increases with h. However, when h — > oo, /ffi (^4 |/t) — > for any finite < a < A, and thus 
i?( T ) (/i, N, A) — > 0; therefore, R^ (h, N, A) should decrease with increasing h when h is sufficiently large. 

C. Rate-Energy Performance Comparison 

As in [4] and [5], there exist rate-energy (R-E) trade-offs in both PS and TS for information and energy transfer. 
R-E trade-offs in PS and TS can be characterized by setting different values of r and A, respectively. Fig. [6] shows 
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R-E trade-offs in PS and TS for N t = 2 and a constant MISO channel h = [1.0 0.56] T For PS, X^ p ) is generated 
by Alamouti code with L = 2 iflQl . For TS, a scalar code cascaded by one single RB is applied when N = 1, 
while the Alamouti code with two RBs is applied when N = 2. The harvested energy normalized by P is denoted 
by E. It is observed that TS yields the best R-E trade-off with N = 1 when E£ < E < h, and with N = 2 when 
E$ < E < E$, while PS yields the best R-E trade-off when < E < E$, where E$ and E$ are shown 
in Fig. [6] Note that at E = 0, i.e., no EH is required as in the conventional MISO system with WIT only, PS 
achieves higher rate than TS since artificial channel fading by random beamforming degrades the AWGN channel 
capacity. However, when the harvested energy exceeds certain thresholds, i.e., Ej? and E$, TS with N = 2 RBs 
and N = 1 RB achieves the best rate performance for a given energy target, respectively. This demonstrates the 
unique usefulness of random beamforming in a multi-antenna SWIPT system even with constant channels. 

It is worth noting that for TS larger information rate is achieved with N = 1 when E$ < E < h, but with 
N = 2 otherwise. This can be explained as follows. For a given h, it can be shown from (fT"3T ) that A — > when 
E ->• h. Thus, with sufficiently small A, we have E^ (h, l,A)m £ (T) (h, 2, A) (note that E^ (h, 2, A) is sightly 
larger than i?( T ) (/i, 1,^4) for small A as discussed for Fig. |4j but the gap between them is negligible as shown in 
Fig. |4] with A = 0.1). On the other hand, with small A, it can be shown from d7J that f^m ( a \h) > Ia\h (° l^ 1 )' 
< a < A, and thus i?( T ) (h, 1, A) > (h, 2, A) from (O, as discussed for Fig. g] Therefore, TS with N = 1 
achieves larger information rate than N = 2 when E is sufficiently large. In contrast, as E — > 0, we have A — > oo 
from (fl"3T i. Then, it can be shown that (h, 1, oo) < R^ (h, 2, oo) since the ergodic capacity of a fading MISO 
channel increases with the number of transmit antennas. Therefore, for TS larger information rate is achieved with 
N = 2 than N = 1 when E is smaller than a certain threshold. 
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III. Performance Analysis in Fading Channel 

In this section, the R-E performances of TS and PS schemes are analyzed in MISO fading channels. It is 
assumed that the MISO channel is independent and identically distributed (i.i.d.) Rayleigh fading channel with 
h ~ CjV(0,Ijv t )» and thus H = ||h|| 2 is a chi-square random variable with 2N t degrees-of-freedom, with the 
following PDF and CDF HI: 

fn (h) = ^h^e-™, (15) 
T(N t ,N t h) 

F n(h) = l- rm . (16) 

In practice, it is possible for Rx to change A for TS or r for PS with the fading MISO channel h for different 
transmission blocks; however, this incurs additional complexity at Rx. In this section, for simplicity, it is assumed 
that A and r are set to be fixed values over different realizations of h. 

A. Achievable Average Information Rate 

We consider that the performance of information transfer is measured by the achievable average rate over fading 
channels. Given N and A, the achievable average rate of TS is denoted by (N,A) = E H [R^ (h,N, A)], 
where R^ (h, N, A) is given by (fl4l) for a given h. However, it is difficult to obtain the closed-form expressions 
for i?( T ) (N,Ays using (0 and CO) for any given N, 1 < N < N t . 

Note that in practice, SWIPT systems usually operate with large transmit power P due to the requirement of 
energy transfer. Therefore, as P is sufficiently large, i?( T ) (h, N, A) in (IT4T > can be approximated as 

(h, N, A) « {A \h) log 2 (P) + Co (h, N, A) , (17) 

where F^ H (A \ h) = j Q A (a\h) da and Co (h, N, A) = j Q A log 2 (a) f^jj (a\h) da, which is a constant not 
related to P. Fig. [7] shows i?( T ) (h, N, A) and its approximation by ([171 ) versus P with N = 1 for different values 
of h and A. It is observed that the approximation in ([171 ) is more accurate as h and/or A increases. It is also 
observed that the gap between the achievable rate and its approximation becomes negligible when P > 30dB even 
with small values of h = 0.1 and A = 0.05. 

With the approximation of R^ (h, N, A) by (TTvT >. we can characterize the asymptotic behavior of R^ (N, A) 
as P becomes large by investigating its pre-log scaling factor, which is given by the following proposition. 
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Fig. 7. Plot of i? (T) (h, N, A) with N = 1, h = 0.1 and 1.0, A = 0.05 and 1.0. 



Proposition 3.1: Given I < N < N t and A > 0, the achievable average rate for TS over the i.i.d. Rayleigh 
fading MISO channel is obtained as (N, A) = (N, A) log 2 (P) + o (log 2 P) with P -)• oo, wherdi 



(iV, A) = lim 



P-+oo log 2 -P 



(18) 



withF^ (a) = E H 
which can be further expressed as 



^A\H ( a I ^ ) denoting the unconditional CDF of A after averaging over the fading distribution, 



AN) 



(o) = l 



^^^ t _ fc (2/3) 



fc=0 



(19) 



where /3 = \JN t Na, and ^ (x) denotes the second-kind modified Bessel function 

7T I-j (a) - jx) 

K 5 {X) = ~ . TZ , 

2 sin (ox) 



with Is (x) denoting the first-kind modified Bessel function 



m=0 



m\T(m + 5 + l) \2J 



X \ 2m+a 



Proof: Please refer to Appendix |B] ■ 
Remark 3.1: In the fading MISO channel, F A denotes the percentage of sub-blocks allocated to ID mode 

for TS. From Proposition 13.11 it is inferred that F A (A) is also the pre-log rate scaling factor of the asymptotic 
achievable average information rate over the MISO fading channel for TS with given A and N. 

Fig. m shows F A N) (A) versus A for TS with Nt = 4 in the i.i.d. Rayleigh fading MISO channel. From Fig. 
[U it is observed that the rate scaling factor F A N ^ (A) for TS increases with decreasing N when A is small, but 



5 / (a;) = o(g (x)) as x — > xo represents that lim 



0, meaning intuitively that / (a;) <C g (x) as x — s> xq. 
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Fig. 8. (A) vs. A when N t = 4. 



decreases with N when A is sufficiently large. As a result, (N, A) scales faster with increasing P for smaller 
value of N when A is small, but scales slower with P when A becomes large. 

On the other hand, the rate scaling factor for PS in the i.i.d. Rayleigh fading MISO channel can be determined 
from © and (031 ) as 

E H [/?( p ) (h.r)] 

(20) 



P->oo log2-P 



r. 



B. Average Harvested Energy 

In this subsection, we study the average harvested energy over the i.i.d. Rayleigh fading MISO channel by TS, 
defined as (N,A) =E H [E^ (h,N,A)], where E^ (h,N,A) is given by O- 

Proposition 3.2: In the i.i.d. Rayleigh fading MISO channel, for given A and N, the average harvested energy 



for TS is given by 



N 



E^ ( N a) =p^-Y 



pN t +k NA 



K Nt . k+l (2/3) , 



(21) 



where f3 and Kg (x) are defined in Proposition 13.11 

Proof: Please refer to Appendix O ■ 
For convenience, we term n( T ) (N, A) = (N, A) jP as the energy scaling factor for TS with increasing 
P. Fig. shows n( T ) (N,A) versus different values of A with Nt = 4. It is observed that the energy scaling factor 
n( T ) (iV, A) for TS behaves in the opposite way of the rate scaling factor (A), as compared to Fig. [8] i.e., 
;q( t ) (JV, A) decreases with N when A is small, but increases with decreasing N when A is sufficiently large. As 
a result, for a given P, E^ (N, A) behaves the same as nm (n,a). 
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Fig. 9. n (T) (N, A) vs. A when N t = 4. 



On the other hand, the energy scaling factor for PS in the i.i.d. Rayleigh fading MISO channel can be easily 
obtained from © and (TT5T ) as 



n<p) 



£ (p) (/i,r) /P=(l-r), 0<r< 1. 



(22) 



The rate and energy scaling factors characterize the asymptotic rate-energy trade-off as P — > oo. Given 1 < 
N < N t , for TS it is easily shown from (fl9l l and (|2T1) that the rate scaling factor A^ T ) (N, 0) = and the energy 
scaling factor U^ T \N, 0) = 1 at A = 0, while A( T ) (N, oo) = 1 and Il( T )(7V, oo) = at A -)■ oo. Note that for TS 
the distribution of the received channel power A (k) at each sub-block becomes different according to N, and as a 
result different asymptotic rate-energy trade-off is achieved when < A^ T ) (iV, ^4) < 1 and < Il( T ) (iV, A) < 1. 
To characterize this trade-off, we have the following theorem. 

Theorem 3.1: In the i.i.d. Rayleigh fading MISO channel, given 1 < N < N t and < A < oo for TS scheme and 
< r < 1 for PS scheme, A( T ) (N, A) > (r) for a given energy scahng factor < nW (N, A) = (r) < 
1; furthermore, given 1 < N < M < N t and < A N , A M < oo for TS schemes, A( T ) (N,A n ) > (M,A M ) 
for a given energy scaling factor < n( T ) (N,A N ) = A( T ) (M,A m ) < 1. 

Proof: Please refer to Appendix |D] ■ 

Fig. [10] shows the rate scaling factor (A( T ) (N, A) for TS and A( p ) (r) for PS) versus energy scaling factor 
(Il( T ) (N, A) for TS and l~[( p ) (r) for PS) with N t = 4. For a given < Il( T ) (N, A) = (r) < 1, the rate 
scaling factor of TS with N = 1, i.e., one single random beam, is the largest among all values of N. In addition, 
A( T ) (N, A) for TS decreases with increasing N, but is always larger than A( p ) (t) for PS. The above observatios 
are in accordance with Theorem 13.11 
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Fig. 10. Rate vs. energy scaling factors with Nt = 4. 



C. Energy Outage Probability 

In this subsection, we study the energy outage probability with a given harvested energy target E at Rx, which is 
defined as PE,out = Pr ( E < E) with E denoting the harvested energy in one block. In particular, we are interested 



in characterizing the asymptotic behavior of PE,out as P — > oo, namely energy diversity order, which is defined as 

, A ,. log p E , out 

d E = - lim — — — - . (23) 

P->oo log-T 

(T) 

Proposition 3.3: In the i.i.d. Rayleigh fading MISO channel, for TS the energy outage probability p Eout with 



P — > oo is given by 



(T) 
PE,out 



, A = 

^N t - (24) 
NA(lnP) ) , A>0. 



Proof: Please refer to Appendix [Gj ■ 
From d23l and (|24l) . it can be verified that the energy diversity order of TS is cf^P = N t when ^4 = 0, i.e., no 
WIT is required, and = with a fixed A > when both WIT and WET are implemented. 



On the other hand, in the i.i.d. Rayleigh fading MISO channels, the energy outage probability of PS with P — > oo 



can be obtained as p^ ut = ( TT^fW ) > ^ — r < ^' t ' lus > from © and the fact that Fjj (h) « /i^* as /i — > 0, we 



obtain the energy diversity order as d^. = Nt, < r < 1. 



It 



Fig. [TT] shows the energy outage probabilities of TS and PS versus P when N t = 2 and E = lOjj It is observed 
that the smallest energy outage probabilities are achieved by TS with A = or equivalently PS with r = 0. When 
A > 0, p E out for TS is observed to decrease slower with increasing P than p E out for PS, since d E = 0, A > 

6 Here, one unit of energy represents the amount of harvested energy in one transmission block when A = 0, P = 1, and h = 1. 
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Fig. 11. Energy outage probability with N t = 2 and E — 10. 

(P) (T) 
for TS while ci^ = Nt, < r < 1, for PS. Furthermore, it is also observed that p E out decreases slower with 

increasing P as A and/or N increases, which is consistent with (|24l) . 
D. Numerical Results 

As inferred from Theorem 13.11 the asymptotic rate scaling factor of TS with N = 1 is the largest among all 
values of N and is also larger than that for PS for a given energy scaling factor as P — > oo. However, it does 
not imply that the largest achievable average rate is always attained for a given average harvested energy when 
P is finite. Therefore, it is necessary to compare the achievable average rates for PS and TS with finite values of 
P. Fig. [12] shows the achievable average rates for PS and TS versus P under the same average energy harvesting 
requirement E (E = E^ (JV, A) /P = £ (p) (r) /P) with JV t = 2. When E = 0.9, the benefit from a larger rate 
scaling factor is clearly observed for TS with N = 1, since it achieves the largest average information rate. When 
E = 0.5, the achievable average rates for TS are similar with N = 1 and 2, but still grow faster with P than that 
for PS. When E = 0.1, the gaps between rate scaling factors of different schemes are small (c.f. Fig. [TOt and as a 
result their achievable average rates become similar. 

Next, Fig. [13] shows the trade-offs between the achievable average rate and energy non-outage probabilities, i.e., 
1 — pE,out, of TS and PS schemes under the same per-block harvested energy requirement E with P = 20dB 
and Nt = 2. When E = 10, it is observed that the energy non-outage (outage) probability of PS is the largest 
(smallest) for a given achievable average rate. However, when E = 50, it is observed that the minimum energy 
outage probability of TS is attained by N = 1 when the achievable average rate is small, but by N = 2 when the 
achievable average rate is larger, while TS with both N = 1 and 2 outperforms PS. 
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Fig. 12. Comparison of the achievable average rate with N t — 2 and the average harvested energy requirement E = 0.1,0.5, and 0.9. 
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Fig. 13. Trade-off between achievable average rate and energy non-outage probability with P = 20dB, Nt = 2, and per-block harvested 
energy requirement E — 10 and 50. 



Remark 3.2: Employing random beamforming at the transmitter requires additional complexity. However, from 
the above results, it is inferred that the achievable average rate is maximized by using only one single RB, i.e., 
N = 1, when P is asymptotically large or at finite P when more harvested energy is required (which is of 
more practical interest). In addition, TS with one single RB also optimizes energy outage performance when P 
is finite and large harvested energy is required in each transmission block. Therefore, TS with one single RB can 
achieve the optimal WET efficiency and/or reliability with a given WIT rate requirement, thus yielding an appealing 
low-complexity implementation for practical systems. 



19 



IV. Alternative Random Beam Designs 



It is worth noting that TS with Gaussian random beams (referred to as TS-G), as considered in the preceding 
sections, may not be practically favorable due to the fact that Gaussian random beams (GRBs) will cause infinite 
transmit power at certain sub-blocks. Instead, artificial channel fading within each transmission block can be 
generated by employing non-Gaussian random beams with constant transmit power for TS. In this section, we 
investigate the performance of TS with two alternative RBs other than GRB, such that the average transmit power 
remains constant within each transmission block, which are given next. 

A. Unitary Random Beams (URBs) 

In this case, N unitary random vectors obtained from the isotropic distribution |[T8l are independently employed 
for the N random beams at the kth sub-block, i.e., n (k), 1 < n < N, \/k. 

With URBs, it is in general difficult to obtain the closed-form expressions for the PDF and CDF of the received 
channel power A (k) at each sub-block conditioned on H = h. However, if we consider the special case of N t = 2 
and N = 1, it is known that with URBs A{k) is uniformly distributed within [0, 2h]. Thus, given a threshold 
A > 0, the amount of harvested energy in each block with URBs can be obtained using (fT2l as 



In the i.i.d. fading MISO channel, the average harvested energy for TS with URBs (referred to as TS-U) given a 
fixed threshold A > is obtained as E^ (A) = J °° E^ (h, A) f H (h) dh, where f H (h) = 4he~ 2h is given by 
(fT5l) for Nt = 2. It is worth noting that in the special case of Nt = 2 and N = 1, the unconditional distribution of 
A (k) with URBs after averaging over the fading channels can be shown to be the exponential distribution, where 




< A < 2h 



A > 2h. 



(25) 



the unconditional PDF is given by (a) = e a . Therefore, E^ (A) can be alternatively obtained as 





(26) 



which is equivalent to E^ (l, 1,^4) for TS-G given by (fT"3T ) with N = 1 and h = 1. Similarly, given A > 0, the 



achievable average transmission rate for TS-U is obtained as 





(27) 
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which is also equivalent to i?( T ) (l, 1, A) for TS-G given by (TT4l ) with h = 1 and N = 1. In addition, given fixed 
A > and per-block energy harvesting requirement i£ > 0, the energy outage probability of TS-U in the i.i.d. 
Rayleigh fading MISO channel with N t = 2 and N = 1 can be obtained from (|25l ) as 

(U ) „ ( e 2 + Ve^ + A^ \ 

Pe'oui = F H I I . (28) 

where Fu (h) is given by ( fT6l ) for JV t = 2. From (1231 and (128T ), it is easily verified that TS-U in the case of Nt = 2 
and N = 1 has the energy diversity order of 0, the same as TS-G. 

B. Binary Random Beams (BRBs) 

In this case, a random subset of N out of Nt transmit antennas at Tx, 1 < N < N t , are selected to transmit at 
each sub-block, which is equivalent to selecting 3> (k) = [4>i{k) 4>2{k) ■ ■ ■ 4 > N {k)] € M. N * xN , \/k, where <j> n (k) = 
[<Pn,i 0) <Pn,2 (fc) • • • (f> n ,N t (k)} T , I < n < N, with n>i (k) G {0, 1}, 1 < i < N t , such that ||0 n (/c)|| 2 = 1 and 
cf>^ n {k)cf) n (k) = 0, n ^ m. We assume that all the subsets of the selected antennas are equally probable. 

Consider the special case of Nt = 2 and N = 1. Denote h = [hi h2\ T , V = max( |/ii| 2 , |/i2| 2 ), and W = 

min( |/ii| 2 , |/i2| 2 )- Note that in this case the received channel power at each sub-block is either A{k) = V or 

A (k) = W, each of which occurs with a probability of 1/2. Thus, given V = v, W = w, and a fixed threshold 

A > 0, the amount of harvested energy in each block with BRBs is obtained using (fT2l ) as 

( P(v + w)/2 , A<w 
(y,w,A) = < Pv/2 , w<A<v (29) 
[ , A>v. 

Similar to TS-U, in the i.i.d. fading MISO channel, it can be shown that with N t = 2 and N = 1 the unconditional 

distribution of A (k) with BRBs after averaging over the fading channels is the exponential distribution, where the 

unconditional PDF is also given by (a) = e~ a . Therefore, given a fixed threshold A > 0, (A) = (A) 

and (A) = (A), where E^ (A) and R^ (A) denote the average harvested energy and achievable 

average information rate for TS with BRBs (referred to as TS-B), respectively, and E^ (Al\ and ^( u ) (yl) for 

TS-U are given by (l26l ) and (l2"7T ). respectively. In addition, given fixed A > and per-block energy harvesting 

requirement E > 0, the energy outage probability of TS-B in the i.i.d. Rayleigh fading MISO channel with N t = 2 

and N = 1 can be obtained from (|29l as (see Appendix |H] for the detailed derivation) 

pgL = (l - e~ A f + 1(A<2D). 2e~ 2 (^) (-1 + e*) (~e A + e 2 ^) 

+1 (A < D) (e~ 2 ( A+D ) (e A - e D ) 2 + e~ A ~ 2D ((-1 + A - D) e A + e D )^j , (30) 
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Fig. 14. Comparison of the achievable average rate for different RB designs with Nt = 2, N = 1, and the average harvested energy 
requirement E = 0.9. 
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Fig. 15. Comparison of energy outage probability for different RB designs with Nt = 2, iV = 1, R = 4 bps/Hz, and E = 50. 



where D = -p, and 1 [x < y) denotes the indicator function given by 



1 (x < y) 



1 , if x < y 



, otherwise. 

From CP, it is shown that both 1 (A < 2D) = and 1 (A < D) = if P > 2E/A, and thus p { ^ out = (1 - e~ A ) 2 . 
Therefore, TS-B in the case of N t = 2 and N = 1 also has the energy diversity order of 0, the same as TS-G and 
TS-U. 

Fig. [14] shows the achievable average rates for TS-G, TS-U, TS-B, and PS versus P for a given average harvested 
energy (normalized by P), E = 0.9. It is observed that the achievable average information rates of TS-U and TS-B 
are the same, which is as expected for the considered case here of Nt = 2 and N = 1. It is also observed that 
the achievable average rates of TS-U and TS-B are larger than that of PS, but smaller than that of TS-G. This 
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result is originated from the fact that the artificial channel fading generated by URBs or BRBs in this case is less 
substantial over time than that by GRBs, due to the limitation of constant average transmit power over sub-blocks 
with URBs or BRBs. 

Fig. Q3] shows the energy outage probabilities of TS-U and TS-B versus P when R = 4 bps/Hz and E = 50, as 
compared to that of PS and TS-G. Among TS schemes, it is observed that the energy outage probability of TS-B 
is the smallest when P is small, i.e., P < 20dB, while that of TS-G is the smallest when P > 20dB. The energy 
outage probability of TS-U is observed to be between that of TS-B and of TS-G. It is also observed that the energy 
outage probability of PS is larger than that of TS schemes when P is small, but is the smallest when P > 22.5dB. 



This paper studied a novel receiver mode switching scheme for the point-to-point MISO SWIPT system when the 
fading channel is only perfectly known at the receiver, but unknown at the transmitter. The proposed scheme exploits 
the benefit of opportunistic energy harvesting over artificial channel fading induced by employing multi-antenna 
random beamforming at the transmitter. By investigating the achievable average information rate, average harvested 
energy/energy outage probability, and their various trade-offs, it is revealed that the proposed scheme yields better 
energy and information transfer performance than the conventional periodic switching without random beamforming 
when the harvested energy requirement is sufficiently large. Particularly, employing one single random beam for 
the proposed scheme is proved to achieve the asymptotically optimal trade-off between the average information rate 
and average harvested energy when transmit power goes to infinity. Moreover, it is shown by simulations that the 
best trade-offs between average information rate and average harvested energy/energy outage probability are also 
achieved by the proposed scheme employing one single random beam for large energy harvesting targets of most 
practical interests, even with finite transmit power. Finally, it is worth noting that the results of this paper can be 
directly applied to the multi-user MISO multicast system for SWIPT, where multiple receivers harvest energy and 
decode (common) information from the signals broadcast by a single transmitter. 



V. Conclusion 



Appendix A 



Derivations of (h, 1,A) and R^ (h, 2, A) 
From with N = 1 and CDS, R {T) {h, I, A) can be expressed as 




(31) 




(32) 
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where (|32l ) is obtained from integrating (1311 by part. By changing a variable as x = 1 + Pa, the integral term in 
d32b can be obtained as 



P 



l + Pa 



-e * do 



'i 

i_ 

e ph 



1 x-l 

—e ph dx 

x 

1 



1 



e ph dx— —e ph dx 

l x ./i+pa x 



1 , *i L ^tt)-*( l inr 



(33) 
(34) 



where, for integer n > 0, E n (z) = e~ zt t~ n dt denotes the exponential integral function; (l34l is obtained by a 
change of variable as y = (l + PA) x for the second integral term in (|33l . From (|3Zj and (|34l , p( T ) (/i, 1, A) is 
obtained as 

«< T ) (», U) - (* (^) - * (i±£)) - e^lo g2 (1 + PI) . (35) 

When N = 2, p( T ) (7t, 2, A) can be derived similarly by integrating (fl4l by part. In this operation, it is necessary 

to apply differentiation of the incomplete Gamma function given by fl9l 

0,0 



^r(N,x)=T(N,x)lnx + xG 3 2 °J ^ 



(36) 



, ^ym,n I a l, ' ' ' ) a P 

where Gr 



T P,Q 



z denotes the Meijer-G function, defined as [Q/7] 9.301] 



(-im,n 

P;9 



ai, • • • , a 

&1, ••• ) 



, , nrfo-s) nr(i-a fe + s ) 

1 / 3=1 fe=l 

27T7 It q P 

*™Jl n r(1 _ 6 . + s) n r(a fc - S 



-z s ds, 



(37) 



j=m+l fc=n+l 

with J L denoting the Barres integral. By the definition of the Meijer-G function, the last term in (l36l can be 



represented by 



xG 



3,0 
2,3 



0,0 



N -1,-1,-1 



1 f T(N- l-a)r(-l-a)r(-l-a) 



2ni 



r - s r - s ) 



x s+1 ds 



G 



3,0 
2,3 



1,1 

0,0, AT 



(38) 



(39) 



where (|39l is achieved by a change of variable as t = s + 1 in (l38l) . By applying (I36l)-d39l to the integration of 
<® by part, (h, 2, A) can be obtained as 



RW(h,2,A) = (^-<r 



( 2(l + PA)\\ 



2 2 / / 2 (1 + PA) \ / 2 



1 J-lrSfil M 
ln2 I 2 ' 3 I o,0,2 



Ph 



G 



3.0 
2.3 



1,1 

0,0,2 



2 (1 + Py4)' 



P/t 



(40) 
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Appendix B 
Proof of Proposition 13. II 

With P — > oo, the achievable average information rate for TS is expressed from (fTTT ) as 



E 



R<W (h, N, A) = E H (A\h) log 2 P + E H [C (h, N,A)]. 



Since Co (h, N, A) is a constant not related to P, so is E# [Co (h, N, A)] and thus E# [Co (/t, iV, A)] is regarded 
as o(log 2 P). 

For an integer N > 1, note that T (JV, x) is equivalently expressed as |[T5l 



N-l 



r(N,x) = (N-l)\e-*J2 
From ©, (H2), and gT), (a) = E H [i^S ("I 71 ) 



,r 



m=0 

is obtained as 



ml 



(a) 



F^ H {a\h)f H (h) dh 



1 — e h 



N-l 

£ 

m=0 



1 fNa\ m \ iVf* 



,N t -l-N t h 



m\\h J T (N t ) 



dh 



Nt N-l 



r (N 



-r y -^(Na) m 
t ) ' ml J 
m=0 0_ 



h N >- m - l e- N > h -^dh. 



(41) 



(42) 



By applying ifTTl 3.471-9] to a, we can obtain ( fT9l . This completes the proof of Proposition 13.11 

Appendix C 
Proof of Proposition 13 .21 

From (Q2), (O, and 0D, # {T) {N,A) = E h [E (t 1 (h, N,A)] is further obtained as 



£ (T) (N, A) = PE h 



N i 

NA x A 1 



PE 



H 



k=0 
N 



k=0 



k\ 



(43) 



By applying IfTTl 3.471-9] to /3 in (l43l . we obtain (T2T1) . This completes the proof of Proposition 13.21 
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Appendix D 
Proof of Theorem |3.1| 

First, the former part of Theorem 13.11 can be proved using the following lemma by considering an arbitrary 
distribution of A with PDF and CDF denoted by gA (a) and Ga (a), respectively, where Ga (a) > for a > 0. 
Lemma D.l: Given Ga (A) = r with < r < 1 and < A < oo, 

POO 

/ agA (a)da> (1 - r) b. (44) 
J A 

Proof: Please refer to Appendix |E] ■ 
For TS, we have the rate scaling factor A^ T ) (iV, A) = F A N ^ (A) from Proposition 13.11 and it can be shown 
from (QUI that F A N) (a) > for a > 0. In addition, the energy scaling factor for TS can be alternatively expressed 



as ll( T ) (N, A) = J A a f A > {a)da with f\ N) (a) = E H (a\h) denoting the unconditional PDF of A after 

averaging over the fading distribution. Furthermore, it can be easily verified that J °° af A N ^ (a)da = 1. Given 
F A N ^ (A) = t with < A < oo and < r < 1, it can thus be verified from Lemma ID. 11 that A( T ) (N, A) + 
n( T ) (N,A) > 1, by substituting b, gA (a), and Ga (a) in Lemma D.l I by 1, f A (a), and F A N ^ (a), respectively. 
Since we have A( p ) (t) + n( p ) (r) = 1 for PS from <gG) and (|22), it follows that A( T ) (N, A) + n( T ) (N, A) > 
A( p ) (r) +n( p ) (r). Therefore, we have A^ (N,A) > (t) for given < n( T ) (N,A) = (r) < 1. This 
proves the former part of Theorem B.lj 

It is worth remarking that Lemma ID. 11 implies that TS in general yields better trade-off between the rate and 
energy scaling factors than PS provided that the average received channel power for TS is the same as that for PS, 
based on which the former part of Theorem 13. II for the i.i.d. Rayleigh fading MISO channel with a fixed threshold 
A is proved. As another example, even for a transmission block with H = h, TS with N RBs, 1 < N < N t , yields 
better trade-off between the rate and energy scaling factors than PS. This can be proved by substituting b, gA (a), 
and Ga (a) in Lemma ID. II by h, f^m {A \ h) in (|7]), and F A ^ (A\h) in (H), respectively. This is originated from 
the fact that for both TS and PS schemes the rate scaling factor is determined by the percentage of sub-blocks 
allocated to ID mode, whereas the energy scaling factor is determined by the percentage of sub-blocks assigned 
to EH mode as well as their channel power values. Note that TS scheme assigns a subset of sub-blocks with the 
largest channel power to EH mode, as inferred from d6). Therefore, given a percentage of sub-blocks allocated to 
EH mode, i.e., 1 — Ga (A) = 1 — r for TS and PS schemes, respectively, the energy scaling factor of TS is larger 
than that of PS while rate scaling factors are the same for both schemes, i.e., Ga (A) = r. 

Next, to prove the latter part of Theorem 13.11 we consider two arbitrary distributions of A with PDFs denoted 
by gA (a) and ua {a), and the corresponding CDFs denoted by Ga {a) and Ua (a), respectively. It is assumed that 
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Jo°° &9A (a) da = J* °° cm a (a) da = b > 0. It is further assumed that Ga (a) > and Ca (a) > for a > 0, and 
Ga (a) and Ua (a) intersect at a = A, satisfying 

G A (a) > U A (a) , if < a < A 

G A (a) = C/a 0) , if a = i (45) 
Ga (a) < ^a (a) , if a > 1. 
Lemma D.2: Given < Ga (A g ) = Ua (A u ) < 1 with < A g ,A u < oo, 

roo roo 

/ agA (a) da > auA (a) da. (46) 
iA g ii„ 

Proof: Please refer to Appendix 10 ■ 
The latter part of Theorem 13.11 can be proved using Lemma ID.2I as follows. Given 1 < N < M < N t for 
TS, it can be verified that f °° a/^ (a) da = f °° af A M ^ (a) da = 1. Furthermore, it can be shown from (|T9) 
that F^ (a) and F^ (a) correspond to Ga (a) and ?7a (a) in d45l ), respectively (c.f. Fig. [8J. By substituting 
/a^ ( a )> /a ( a )' ^A ( a )' an d F^-* (a) for $a («), «A (a), Ga (a), and Ua (a) in Lemma |DTl respectively, it 
can be verified that n( T ) (N,A N ) > (M,A M ) when A( T ) (N,A n ) = (M, Am), < A N ,A M < oo 
since A( T ) (N,A n ) = F^ V) (a) from Proposition ED and IlW (N,A N ) = Jj^ af A N) (a)cfo. This guarantees that 
A( T ) (N,A n ) > AW (M,A m ) for given < n( T ) (N,A N ) = (M,A M ) < 1, since both A< T ) (N,A N ) 
and A^ T ) (M, Am) decrease monotonically with increasing n( T ) (N, An) and I[( T ) (M, Am), respectively. This 
proves the latter part of Theorem 13.11 

As a remark, Lemma ID.2I compares the trade-offs between the rate and energy scaling factors in TS schemes 
with two different distributions of channel power induced by different values of N provided that both distributions 
have the same average channel power and satisfy the condition in d45l ). The latter part of Theorem 13. II for the i.i.d 
Rayleigh fading MISO channel with a fixed A is one application of Lemma ID.2I As another example, even for 
a transmission block with H = h, better trade-off between the rate and energy scaling factors is attained with N 
than M RBs, 1 < N < M < N t , since F^ (a \h) and F^ (a\h) correspond to G A (a) and U A (a) in d45l 
respectively, as shown from ©. This is due to the fact that the artificial channel fading is more substantial when 
smaller number of RBs is employed, and the argument similarly as for Lemma |D.1| 
Combining the proofs for both the above two parts, Theorem 13. II is proved. 

Appendix E 
Proof of Lemma IdTTI 

Integrating by part, agA (a) da can be evaluated as 

i-A pA 

\ ag A (a) da = AG a - G A (a) da. (47) 
Jo Jo 
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Assume that A is given such that Ga (A) = r, < r < 1. From (|44T > and (14"7T ), we have 

/■oo 

E = agA (a) da — (1 — r) 6 
./A 

/■oo /-A 

= / ag A (a)da- / a&4 (a) da - 6 (l - G A (A)) 
Jo Jo 

= {b-A) G A {A) + / G A (A) da, (48) 
Jo 



and the resulting derivative 



, = ^ = (b-A)g A (A), (49) 



where \i = is achieved at A = b. 

From (l48l and (l49l . it is observed that i? > with < A < b, and increases monotonically with A until 
A = b. For A > b, it is observed that /i < and thus E decreases monotonically with increasing A. Because 
lim f?° aoA (a) cia = 0, it can be shown that E > even with A > b. Therefore, it is verified that E > in 

A— >oo 

for any ^4 > 0. This completes the proof of Lemma ID. II 

Appendix F 
Proof of Lemma ID~2l 



Denote A = Ga(A 9 ) = Ua(A u ), < A g ,A u < oo. From (g6j and (l47i we have 

E = agA (a) da — auA (a) da 
J A g JA U 

= / auA (a) da — agA (a) da (50) 
Jo Jo 

rA 9 pA n 
= A(A U -A g )+ G A {a)da- U A (a)da. (51) 
Jo Jo 

According to (I43T ). there are three cases addressed as follows for given < A < 1. 

1) A g = A u = A: According to d45i Ga(A) = Ua(A) = A. Since it is assumed in that Ga(cl) > U A (a) 
with < a < A, E is evaluated from (IBTI ) as 



E = [ (G A (a) - U A (a)) da > 0. (52) 
J o 

2) < A g , A u < A: It can be inferred from (@3) that ^ 9 < A u < A, which results in < A < A. From (ED, 
we have 

£ = U A (A u ) (A u - A g ) - U A (a) da + (G A (a) - U A (a)) da . (53) 

J A„ Jo 



2S 



Since A g < A u and Ga{o) > Ua{o) with < a < A, it can be verified that a > and /3 > 0, and thus E > 0. 

3 J A < A g ,A u < oo: It can be inferred from (@3]l that 4 < A u < A g , which results in A < A < 1. From (|50l >. 
we have 

E 1 = J a (ua (a) — gA (d)) da — (j agA{a)da — j' auA{a)da^j, (54) 

V V ' V v ' 

= 8 A 

= £ 

with 5 > as shown in (I50l)-d52l). In addition, it can be verified that 

lim e = lim / clqa (a) da — Jim / auA (a) da = 0, (55) 

A^A Ag-^AJA A u -^AJa 

i- A i- A 

lim e = lim e= / a-u^ (a) da — / agA{a)da = 5. (56) 

A u — >oo 

Since zfe = dfc^ (4) = <M (4) and £ = gj-t^ (A) = tm (A), we have 

d 1 d r*» , , . i d ^ 



-e 



ag A (a)da- dA / aw^ (a) da 



dA di; si: dAu 

= A g - A u > 0. (57) 

From (l54l-(l57Ti. it can be verified that E monotonically decreases from 5 to with increasing A with A < A < 1, 
i.e., E > with A < A g , A u < oo. 

Combining the above three cases, Lemma ID.2I is thus proved. 

Appendix G 
Proof of Proposition [33] 

Given a transmission block with H = h, E^ (h, N, 0) = Ph from ([13]). In the i.i.d. Rayleigh fading MISO 
channel, given E > 0, p ( 5> for TS with A = can be approximated as P — > oo by lim Pr I h < Jj I = 

' P— »oo V / 

lim F/f ( 4 ) = ( i ) \ since lim F# (h) = h Nt in ([16]). This proves the first equality in (|24) for A = 0. 

When A > 0, the harvested energy per block for a given h, E^ (h, TV, A) in (QjJ, is a monotonically increasing 
function of h, since T (a, x) is a monotonically decreasing function of x. For a given energy requirement E > 0, 
denote a as the minimum value of h such that i£( T ) (h, TV, A) > E, i.e., 

r(jv + i) "p- (58) 

Thus, the energy outage probability for TS with given TV, A, and 22 is obtained as pffi t (N,A,E) = Fh(Ji). 
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Since 6 (h) increases with h, with P — > oo it then follows from (1581 ) that h — > 0, under which we have 

»(s)-k^>£i(^y (59) 



fc=0 

where (|59l l is obtained from d4ll since A^ > 1 is an integer. Since ln# (h) = In from (l60l ). we have 

/ FAM \ 

(JV-l)lnx- NAx = ]xl\ w , (61) 

where x = 1/h. With /i — >• 0, i.e., x — > oo, the left-hand side of (f6Tb can be further approximated as — NAx. 
Therefore, h as P — > oo can be approximated by 

h = NA\ In | \ ' J) « iV^(lnP) -1 . (62) 

From (l62l and the fact that F# (/i) /i^* as h — > 0, we obtain the second equality in d24l) for A > 0. 
From the proofs for both the first and second equalities in (|24T) . Proposition 13.31 is thus proved. 

Appendix H 
Derivation of ([301) 



From (|29l l, the energy outage probability of TS-B is obtained as 

PeIui = Pr i v < A ) + Pr y» < A < v > ^ < E) + Pr ^ > A, P W) < eJ . (63) 
First, Pr [y < A) in (l63l is given by 

Pr (i; < .A) = Fy (A) = (l - e~ A ) 2 , (64) 



where Fy (v) denotes the CDF of V = max( |/ii| 2 , |/i2| 2 ), given by Fy (v) = {1 — e v ) , since both \h\\ 2 and 

2 

\h2\ are independent exponential random variables. 
The second term in (l63l can be obtained as 



Pr fw < A < v, ^ < F?j = Pr (w < A, A < v < 2D, w < v) 

I-2D rA 

= l(A<2D) / f v ,w (v, w) dwdv 

J A JO 

= 1(A< 2D) 2e~ 2 ( A+D ) (-1 + e A ) (-e A + e 2B ) , (65) 
where D = p and fv,w (v, w) denotes the joint PDF for V and W given by fy,w (v, w) = 2e~ v e~ w , v > w. 
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Finally, the last term in (|63l l can be obtained as 

/ , P(v + w) A 
Pr lw > A, v ; < E\ 

= Pr (w > A, v + w < 2D) 

(pB rv p2B-A j-2B-v \ 

/_ /_ fv,w (v, w) dwdv + / / /y;w («, u>) dwdv 

J a J A Jb J a J 

= 1{A<D) (e- 2 ( A+D ) (e A -e D )\ e" A " 2D ((-1 + A - D) e A + e D )) . (66) 
From (l63t-(l66l. we can obtain (|30T) . 
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